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Abstract 

We investigate particle production near extra species loci (ESL) in a higher dimensional field space and 
derive a speed limit in moduli space at weak coupling. This terminal velocity is set by the characteristic 
ESDseparation and the coupling of the extra degrees of freedom to the moduli, but it is independent of the 
moduli's potential if the dimensionality of the field space is considerably larger than the dimensionality of 
the loci, D 3> d. Once the terminal velocity is approached, particles are produced at a plethora of nearby 
ESLs, preventing a further increase in speed via their backreaction. It is possible to drive inflation at the 
terminal velocity, providing a generalization of trapped inflation with attractive features: we find that more 
than sixty e-folds of inflation for sub-Planckian excursions in field space are possible if ESLs are ubiquitous, 
without fine tuning of initial conditions and less tuned potentials. We construct a simple, observationally 
viable model with a slightly red scalar power-spectrum and suppressed gravitational waves; we comment 
on the presence of additional observational signatures originating from IR-cascading and individual massive 
particles. We also show that moduli-trapping at an ESL is suppressed for D d, hindering dynamical 
selection of high-symmetry vacua on the landscape based on this mechanism. 
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I. INTRODUCTION 



The presence of many light fields, or moduli, is a common feature of string theory. At late times, 
the expectation values of these fields determine low energy observables; hence, their evolution is 
heavily constraint from the time of nucleosynthesis, although they are expected to be dynamical in 
the very early universe. For these reasons, moduli trapping is an important aspect of inflationary 



model building in string theory, see 



4_5j] for reviews; often all but one degree of freedom are 
stabilized by construction during inflation, but how they stabilize is seldom addressed. 

A possible dynamical stabilization process is the String Higgs effect P, Q], which singles out 
certain locations in moduli space. At these loci additional degrees of freedom become light, often 
due to the presence of an enhanced symmetry, and are produced ?| if the moduli approach 
the location. The process of particle production is identical to the one examined in preheating 
see 

for reviews. Backreaction of these new states on the moduli can be dramatic: 



m „duli can betrapped flflQW or, drive inflation, th. inflaton velocity can decrease 



temporarily [l6|], as in trapped inflation 



can be found in a system of parallel D-branes 



lal. A concrete realization of the string Higgs effect 



191 ]. whose separation is a modulus field from the 



four-dimensional point of view. If they come close to each other, strings stretching between the 
branes become massless, gauge symmetries are enhanced and the branes stick to each other due to 
backreaction - the modulus is trapped. For a small selection of other examples, see 20-4271] . 



In this paper, we examine the consequences of particle production near such extra species loci 
(ESL) if the dimensionality of moduli space D is large compared to the dimensionality of the 
locus, d. D > 1 is natural in string theory, leading to the notion of a landscape [28|, yet the 
effect of ESLs has primarily been investigated for D — d = 1,2 [?], Q, [29]. Based on geometric 
arguments, we show in Sec. [II] that trapping is suppressed if the dimensionality of moduli space is 
larger than the dimensionality of a given locus (D > d + 1). This result is expected, since there is 
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no classical attraction towards ESLs and it is improbable to run head on into an ESL if D 3> d. 
This means that the mere presence of ESLs does not guarantee a dynamical preference of high 
symmetry states for moduli 



30, 



311 ] . However, in Sec. lIIIB II we show that the presence of many 



loci with a characteristic inter-ESL distance x leads to a gener al sp eed limit, or terminal velocity, 



321 ] . leading to DBI-inflation 



33|. 



on moduli-space. At strong coupling, speed limits are known 
Here, we derive a speed limit at weak coupling caused by the combined backreaction of particles 
produced at many ESLs in the vicinity of the trajectory. 

We take a bottom up approach, assuming the viability of low-energy effective field theory 
and treating the characteristic distance between ESLs, as well as their dimensionality, as free 
parameters. Furthermore, we model the additional light degrees of freedom by a massless scalar 
field that couples to the moduli via interactions of the type g 2 x 2( f 2 i as in [8j]. In this notation, the 
speed limit takes the simple form |<^| < gx 2 = vt, given that the classical trajectory is reasonably 
straight (we allow for a classical potential for the moduli)and D is large. The presence of a terminal 
velocity offers the opportunity to drive inflation by moduli with hitherto unsuitable potentials. 

This type of inflation is a generalization of trapped inflation 0, which recently re-surfaced 
as monodromy inflation [l^ ]. but with several crucial differences: firstly, particle production at any 
given ESL is minor because ESLs are not approached too closely, in contrast to the one-dimensional 



case m 



17J; there, extra species points (ESPs) are encountered head on and all particle production 



occurs at a single ESP at any given time. Secondly, in the large D limit the velocity at which 
the trajectory is traversed is independent of the slope of the potential as long as the slope is steep 
enough. Again, this differs from [171 ]. where the slope still determines the speed and observational 
parameters, such as the scalar spectral index. 

We investigate this type of trapped inflation in higher dimensions in Sec. IIIICl where we en- 
counter several attractive features: focusing on a model with quadratic potentials for the moduli, 
more than sixty e-folds of inflation can be achieved with sub-Planckian field excursions if v% is small 
enough. The inflationary scale is set by the COBE normalization, but the dynamics is not sensitive 
to the actual shape of the potential as long as a single, reasonably straight classical trajectory with 
large slope is present. This is, in a sense, opposite to setups needed for slow-roll inflation, where 
the potential needs to be shallow over long ranges in field space. As a consequence, the ^-problem 
is alleviated. Furthermore, no fine-tuning of the initial speed is needed: if the initial speed is 
large, particle production is strong, causing strong backreaction and temporary trapping; after the 
produced particles are diluted away due to the expansion of the universe, the speed picks up again 
and the terminal velocity is approached from below. In this sense, trapped inflation is an attractor 
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solution. 

The ongoing particle production during inflation has several additional observational conse- 
quences, Sec. IIII D 1 1 and IIII D~2l Since x-particles scatter off the inn at on condensate, they cause 
IR-cascading and additional contributions to the power-spectrum 



34 



35]; furthermore, x-particles 



36l438||. We plan to investigate 



dilute and become massive once the trajectory moves away from their ESL of production. Thus, 
one needs to consider the effect of individual massive particles that are coupled to the moduli, 
which could lead to additional circular cold-spots in the CMBR 
both effects, as well as the generation of non-Gaussianities, in a future publication, accompanied 
by a more rigorous string theoretical implementation of trapped inflation in higher dimensions. 

Readers familiar with the notion of particle production, as discussed in [8J, and a primary 
interest in the terminal velocity and consequences for inflation may skip Sec. [TT] and go directly 
to the main part of this work, Sec. IIII1 On the other hand, readers with a primary interest on 
a dynamical selection principle on the landscape will find Sec. [TT] and Appendix [C] useful (these 
sections review Q] and contain minor extensions and consequence which are implied to in jg| but 
not fully spelled out). 

II. QUANTUM TRAPPING IN MODULI SPACE 

When are moduli affected by the presence of extra species loci? After estimating the distance 
that moduli can move classically in Sec. Ill Al we review the notion of quantum moduli trapping in 
Sec. Ill Bl as developed in [8] and estimate the trapping probability in Sec. Ill CI We comment on 
drifting moduli in Sec. Ill Dl and Appendix O Readers familiar with 8j may want to skip Sec. Ill Al 
and Sec. HlBl 

A. Length of the Field- Trajectory 

Let's start by reviewing the effect of Hubble damping on the evolution of moduli fields, as 
discussed in [8]. Consider a D-dimensional moduli space with freely moving fields (<pi, i = 1 . . . D), 
that is, without a potential that could trap the fields classically. Additionally, assume a flat moduli 
space metric, leading to a straight trajectory. Defining (p c i as the effective field along this trajectory, 
its equation of motion is 

(p c l + 3Hip cl = 0, (1) 
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where a dot denotes a derivative with respect to cosmic time and H = a/ a is the Hubble parameter. 
For a constant equation of state parameter w = p/p the Hubble parameter becomes H = /3/t with 
(3 = 2/(3 + 3w), and (P) can be integrated to 

3/3 



(t) > 



where t> = </i c ;(io) > 0. 

If the universe is dominated by freely moving scalar fields, that is if w = 1, the length of the 
field trajectory 



increases logarithmically 



t 

On the other hand, if w < 1 that is (3 > 1/3, equation ([2]) yields 

3/3 



\<Pd(t) - Vci(to)\ (3) 



8 = «*dM^r) . ( 4 ) 



3(3-1 \t ) + 3(3-1 ' ® 

which is bounded from above by s ^ vto/(3(3 — 1). Since w < 1, we have u 2 /2 < 3M^H 2 so that 
s < ,/6PM p /(3P - 1). 

To get intuition for the expected values of s in realistic scenarios, let's consider concrete numbers; 
in an inflating universe with (3 — > oo the upper bound becomes s < \^6M p /3, less than M p regardless 
of the initial velocity. However, if (3 is closer to 1/3, the length of the trajectory can increase above 
M p . For instance, consider a phase dominated by kinetic energy, 

with e < 1, where a small contribution to the total energy density p redshifting as p c (t) = 
ep(to)(to/t) 1+Wc (slower than v 2 oc a -6 oc t~ 2 ) is also present. p c takes over at t c = tQS~ 1 /( 1 ~ Wc ) 
when v 2 (t c )/2 = p(t c )/2. Using ([4]) and ([5]), we get a total path length in the limit t — > oo of 

s = -—^—vtoln(e) + V3M pl — . (7) 

1 - w c ^ 3(1 - w c ) 

Thus, an extended kinetic phase is needed to achieve super-Planckian values for s. For example, 
if moduli have an initial Planckian kinetic energy (v ~ s/2M 2 ) and if they dominate until high- 
scale slow-roll inflation sets in around Hi n j ~ 10~ 5 M p (w c ~ —1), we get e ~ 3 x 10~ 10 and s ~ 
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— -v/l /6M p ln(e) ~ 9M p . Conceivably, moduli could also dominate after inflation ended i.e. down 
to the SUSY breaking scale p e ~ TeU 4 , so that e ~ 3 x 10 61 and s ~ O(10 2 )M p . 

In summation, if the kinetic energy of the moduli does not dominate the energy density of the 
universe, the moduli travel a limited stretch in field space until they come to rest due to Hubble 
friction, s ~ M p ; for instance, they do not travel far during inflation or in a radiation/matter 
dominated universe . If the moduli's kinetic energy dominates, i.e. before a standard high scale 
inflationary phase or in the interval after inflation but before SUSY breaking, they can travel 
farther (s boforo inf . < 10M p and s aftcr inf . < O(10 2 )M p ). In the following sections we keep in mind the 
interval 

1 < s/M p < 10 2 . (8) 



B. Quantum Trapping 

Even in the absence of a classical potential, moduli can be trapped near locations where addi- 
tional particle species become light. These additional degrees of freedom are produced quantum 
mechanically if the field trajectory gets close enough to an extra species' locus (ESL), and their 
backreaction can force the trajectory in an orbit around it. Since Hubble damping in an expanding 



re^origin of the additional 
, so that ESL or ESP 



universe causes an inward spiraling trajectory, the moduli get trapped. T 
light degrees of freedom are often enhanced symmetries at the locus 
can also denote enhanced symmetry locus or point. Quantum moduli trapping was proposed in 
[3, ll| and was recently investigated in 29 ] . A field theoretical description of this phenomenon is 
identical to the resonant particle production after inflation, usually referred to as preheating jg|.[lo|. 



see i.e. [111 . [12j for reviews. Our first goal is to estimate a characteristic distance to the ESL below 
which moduli trapping is common. 

Consider a single ESP at the origin where additional light states appear. We model this situation 
by introducing an new scalar field \ coupled to the freely moving moduli via a quadratic interaction 

L = \ E d^id^i + \d^X ~ y X 2 E > (9) 

i=l i=l 

without any bare x mass, as in [8]. Classically, x = along with a straight trajectory in moduli 
space is a solution. Starting with this solution, the trajectory and the ESP span a plane in moduli 



1 They could also be displaced by quantum mechanical fluctuations during inflation, see Appendix IC 1 
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No trapping (io>|i Trapping |io<|J, 
?(t) 





FIG. 1: Schematic of a trapping event: if the trajectory <p(i) comes close to an ESL (gray circle), that is if 
the impact parameter /iq is smaller than the ESLs effective radius /i < y/v/g, light particles are produced; 
their backreaction on ipi, along with Hubble damping, can trap the trajectory in the vicinity of the ESL. 



space. Hence, we can redefine the fields so that the straight trajectory takes the simple form 

tpi = vt (10) 

¥>2 = MO (11) 

ifi = for i > 3 (12) 

where we ignored Hubble friction so that v = const. Hence, to discuss the effect of a single ESL, 
the problem reduces to the two dimensional case investigated in [8]. Efficient particle production 
in a Fourier mode Xk occurs if the non-adiabaticity parameter 



D 



1/2 



(13) 



satisfies |1 



oj- 



8=1 



> 1. 



(14) 



Since the effective mass of the x~ ne ld is given by rntit) = g 2 Y2i V 3 ? (*) we see that this requirement 
is satisfied if the x-fi e ld is light. Once x _ P ar ticles are produced, they induce a classical confin- 
ing potential for the moduli, which gets stronger with successive bursts of particle production. 
Evaluating (|14j) at the point where the ESP is closest, we conclude that trajectories with impact 
parameters smaller than 



Mo ^ 



= u 



(15) 



have a strong chance to get trapped near the ESP (see Fig. [T] for a schematic). The value of g, 
and correspondingly the value of the characteristic impact parameter fx below which trapping is 
common, is model dependent. In the case of a higher dimensional ESL the discussion is analogous, 
with no denoting the minimal distance between the straight classical trajectory and the ESL. Since 
fx oc y/v, the value of the characteristic impact parameter decreases drastically once fields slow 
down. 

The actual trapping event is complicated by Hubble damping as well as the fragmentation 
of fields once backreaction becomes important. Except in the simplest one-dimensional cases, a 
proper discussion of trapping requires numerical simulations, as in preheating 3914431] . 



In this section, we take a binary point of view and assume that the moduli get trapped once 
they approach an ESL to within a distance smaller than fx (we go beyond this binary treatment in 
Sec. IIIII) . We further treat fx as a constant 2 . 



C. Trapping Probability 

As we saw in section Hi A\ Hubble damping limits the overall path length to s ^ few x M p i in 
an expanding universe. If the field trajectory <p(t) approaches an ESL to within the characteristic 
impact parameter fx, trapping of the moduli is common, see Sec. IIIBI and Fig. [TJ We denote the 
dimensionality of the ESL by d (d < D), so that d = corresponds to an ESP, d = 1 to a line, etc. 

Further, we assume that ESLs are spread over moduli space with a characteristic distance x (we 
make this more precise below). For x > fx, the probability of trapping is less than one, provided that 
the <fi start out at a random initial position and velocity. Our goal is to estimate the dependence 
of this probability on the ratios s/x and fx/x, as well as the scaling with the moduli space's and 
ESLs' dimensions, D and d. Throughout this section, we take s > x > fx (see Sec. IIIIBl for x ~ fx). 



1. Extra Species Points (ESP, d = 0) 

We assume that ESPs are common and distributed over field space with a number density of 
a = l/x D . Consider a D-dimensional sphere of radius s centered at <p(0), enclosing the region that 



2 We will find that ptrap oc [iP~ d ~ 1 ; hence, trapping becomes less likely as v decreases due to Hubble friction, so that 
we overestimate ptrap by treating /i as a constant over the whole trajectory. A better estimate could be achieved 
by integrating a probability density over the trajectory (p(t) using u(t). Since the latter is model dependent, we 
will be satisfied with the overestimate caused by treating u ~ const (the correction usually introduces factors of 
order one). 
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D =2 dimensional Moduli Space 



♦i 



FIG. 2: Schematic of a higher dimensional moduli space (D=2 is shown): the large light- gray sphere of 
radius s encloses the accessible volume by <p(t). The small dark-gray spheres are ESPs (d = 0) of effective 
radius /i, distributed with an average inter-ESP distance x. 



the moduli can access, see Fig. [2j In this sphere we find 

TrD Vds D 

n-ESL = <rV D s = -p 

ESPs that are potentially reachable by <p(t). Here 

„D/2 



v D 



(16) 



(17) 



r(f + i) 

is the volume of a D-dimensional unit-sphere, and T is the Gamma-function (r(x) ~ 
\/2ttx x ~ 1 / 2 exp(— x) for large x). Since the majority of the volume in a sphere is close to its 
surface for large D, the distance from ty?(Q) to the majority of the accessible ESPs is given by a 
distance close to s, the radius of the sphere. Thus, we approximate the distance of all ESPs to the 
center of the sphere by s. A lower bound for the probability of trapping can then be estimated by 
the ratio of the sphere's surface blocked by the ESPs to the total surface area 

nESpVD-iH ^ 1 



Pel ~ 



Sd-is - 1 



where 



5 D _! = DV D 



(18) 



(19) 



is the surface area of a D-dimensional unit-sphere, see Fig. [3] for a schematic. Here, we assume 
that ESPs are sparse enough to prevent overlapping of blocked regions on the sphere's surface so 
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that our results are only valid for p c \ -C 1. Further, we take D > 2, since the one dimensional case 
is trivial: the probability of getting trapped is p c [ = 1 if x < s and p c i = s/x for x > s. (TTSj) is a 
lower bound, since even a few ESPs close to the center of the sphere, which we treated as if they 
where sitting on the surface, can block out a large fraction of the surface. We come back to this 
issue in appendix 1X1 

With (|16p the lower bound on the probability reads 

> 1 T7 s f^ - 1 



»« ~ D^z If) ' < 2 °» 



In line with expectations based on [8], the more ESPs are within reach, the higher the probability 
of trapping. However, since /j, < x and p c i oc (fi/x) ^ 1 we see that trapping becomes heavily 
suppressed for D > 2; thus, for large D, the starting point <^(0) needs to be within close range of 
an ESP, that is within a distance of order fi or less, to have a reasonable chance to get trapped. 
This should be compared to the one dimensional case, where the relevant distance is s/x: the 
modulus field is likely to be trapped as long as a single ESP is within reach. 

The suppression p c \ oc (fi/x)®' 1 has important consequences for vacuum selection on the land- 
scape. When the notion of quantum moduli trapping near ESPs was investigated in ?[ a 
dynamical selection principle to single out certain vacua associated with additional symmetries on 
the landscape, as used in 



30, 



311 ]. was a strong motivation. However, for large D we see that (p 
needs to start out or be accurately aimed at an ESP in order to get trapped near it, hindering this 
dynamical selection. This argument remains true for loci of higher dimensions as long as D S> d 
and x > [i. 



2. Extra Species Loci (ESL, general d) 

There are many cases where symmetries are enhanced (and additional light states appear) not 
at points, but at loci of higher dimensionality. We can easily generalize the derivation of the 
last section to arbitrary d. We first assume a distribution of the ESLs on moduli-space with a 
characteristic distance x, similar to the case of ESPs; for example, if d = 1, one can visualize the 
loci as a string network with characteristic inter-string distance x. Thus, we can define an analog 
to the number density in the ESP case via a = l/x D , so that the blocked (by ESLs) surface area of 
the D-dimensional sphere with radius s is bigger than Vr>s D ax d /j> D ~ 1 ~ d V£)-i-d- Here we treated 
all loci as if they were located at the surface of the sphere (hence, we find again a lower bound 
for the probability) and that they are sparse enough to ignore overlapping. As a consequence, the 
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\r M D-l-d d 

p V °-'-^ x V D s D o 



FIG. 3: Visualization of the trapping probability (the case D — 2 and d — is depicted): we first push all 
ESLs (dark-gray) towards the boundary of the light-gray sphere of radius s. The ratio of the area that is 
blocked by the ESLs to the total surface area gives a lower bound of the trapping probability. 



trapping probability satisfies 



> 1 T/ S 

Pel <j -prKD-l-d- ( - 



D-X-d 



(21) 



For d < D — 1 we observe a suppression of the trapping probability oc (fi/x) D ~ 1 ~ d . Furthermore, 
Vo-i-d ->0asD-l-d-> oo, providing yet another suppression. See Appendix [B] to see how 
this result applies to two simple examples. 

The factor of 1/D is absent if we estimate the probability more accurately by means of a volume 
argument, see Appendix [A] which does not require any artificial shifting of ESLs, yielding 



l-d- 



S (11 



X \x 



D-l-d 



(22) 



Here, iVj_ depends on the distributions of ESLs; for example, the grid-like distribution in Appendix 
|A] leads to (El, 



N ± 



D - 1 
d 



(D-l) 



d\{D-l-d)\ 



(23) 



D. Comments on Drifting Moduli 



During inflation any classical movement of moduli is quickly damped by Hubble friction, yet 
quantum mechanical fluctuations can still cause fields to drift through moduli space given that 
Hubble induced corrections to the moduli's masses are negligible and no classical potential is 
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present. If we follow the field value in a given Hubble patch, we can observe a random walk. 
Hence, one might expect an increased trapping probability, since a random walk is more volume 
filling that a straight trajectory. Since the Hausdorff dimension of a random walk is 2, that is 
a random walk is area filling, we expect the trapping probability to be slightly less suppressed, 
p rw oc ([i/x) D ~ 2 ~ d instead of p c i oc (fi/x) D ~ 1 ~ d . This expectation is correct 3 , see appendix IC 31 
We note that if the fields perform a random walk as described in appendix [UJ and trapping does 
occur, then the trapping event has to take place before the last sixty e-folds of inflation; in that 
case, our entire Hubble patch today lies in a region that was inside the Hubble patch at the time 
of trapping and as a result, the moduli fields have the same value everywhere in our observable 
universe. 

However, we are primarily interested in the limit D 3> d, leading to a strongly suppressed 
trapping probability if \i < x. Then the mere presence of ESLs does not guarantee late time values 
of moduli fields near these locations; see appendix IC 41 for more details. 

In section fill Bl we investigate the case \x ~ x, in which case we expect particle production to 
be important; we further allow for a classical potential for the (fi, so that they may be responsible 
for inflation. 



III. MODULI AS INFLATONS 



So far, we focused on fields whose potential energy is negligible compared to both, their kinetic 
energy and 3H 2 /(8itG) (from here on we set 8ttG = 1). If the potential energy dominates, a 
phase of moduli-driven inflation results (see 5[ for a recent review on inflationary models in string 
theory); successful moduli trapping goes hand in hand with (p)reheating, while failed trapping 
events can still have interesting consequences. 

If moduli serve as inflatons, a successful trapping event is unwanted if it occurs during the last 
sixty e-folds of inflation. However, failed trapping events during inflation, that is some particle 
production at an ESL, slow down the fields and render steeper potential viable for slow-roll inflation. 
This is the idea behind trapped inflation |s], [jjj], recently resurfaced as a special case of single 



field monodromy inflation 



18] with sub-Planckian field excursions. Here, ESPs are repeatedly 



encountered by the inflaton as the angular modulus in a string compactification revolves around 



3 Note that /i in appendix IC 31 has a different meaning: it is not related to the speed of the fields, but it is set by 
the distance at which the extra degrees of freedom can be described by a massless scalar in field theory. However, 
the scaling p rw oc (fi/x) D ~ 2 ~ d is also valid for fi — \fvjg if a classical trajectory is intertwined on small scales so 
that it can be approximated by a random walk in D dimensions. 
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a periodic direction 



18j], see also 4^]. Due to the classical (smooth) potential, the inflaton does 
not get trapped for long (even though D = 1 and d = 0) but it slows down until inflation dilutes 
the produced particles sufficiently and the field picks up speed again. After leaving the ESP, the 
produced ^-particles quickly become heavy, re-scatter off the homogeneous inflaton condensate and 
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34, 



These fluctuations can dominate over the 



351 ] . such as a bump in the power-spectrum 



cause a cascade of inflaton fluctuations in the IR 
primordial ones and lead to observable signatures 
for each failed trapping event and potentially large non-Gaussianities. 

One may wonder if it is conceivable to increase the number of fields in trapped inflation, for as 
long as each of the inflatons comes close to ESLs during inflation, the above mentioned effects will 
be present. 



A. x fi, Suppressed Trapping 

Since the probability of getting close enough to ESLs depends on the type of movement through 
moduli space, we need to discern between three classes of inflation 

1. Fields with a smooth potential. A smooth classical trajectory in field space results, without 
any tunneling events. 

2. An overall smooth potential with small random bumps and dips. The classical trajectory 
can be described by a biased random walk; the overall drift is determined by the averaged 
potential, while the superimposed Brownian motion is caused by scattering events on the 
bumps, see i.e. 

3. A potential with prominent valleys and mountains. The evolution of fields is the result of 
repeated tunneling/scattering events, not slow-roll, see i.e. j^f], Q|- 

If D > d + 1 and fields roll slowly, as in case one, we expect trapping events to be rare, since 
Pd oc {n/x) D ~ 1 ~ d from (|52]1 and \i = y/v/g from (fl~5j) . 

If the trajectory is intertwined, as in case two, we need to discern between the drift velocity v, 
which is again small, and the velocity along the trajectory v = (^ t£>|) . These two velocities 
are comparable if small scattering events dominate; thus, we expect again an absence of moduli 
trapping if D > d + 1. However, if v 3> v, for example due to frequent large angle scattering 
events, we note that ji is not suppressed. The trapping probability of the resulting biased random 
walk lies between the classical one in (I22p and p rw oc (fi/x) D - 2 - d from (fC27l) . To be concrete, if 
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vN/H 3> sV Nn, where s is the mean free path between two scattering events (not H/2tt), N is 
the number of e-folds and n is the number of scattering events per e-fold, we may approximate 
Pbrw ~ Pel- On the other hand, if the classical drift is negligible, vN/H -C V Nns, we may use 

Pbrw ~ Prw 

In case three, the trajectory is well described by an unbiased random walk so that (|C27|) is 
applicable. However, the step width s, number of steps per e-fold n, and fj, need to be reconsidered, 
because the trajectory is the result of repeated tunneling/scattering events. 

Regarding the viability of trapped inflation in the presence of multiple fields, we note that if p c i 
is applicable, there needs to be a series of closed D—l dimensional ESL-hypersurfaces encompassing 
the origin to which <p(t) is heading; since d = D — l, the trapping probability ([22]) is not suppressed, 
particle production can occur and trapped inflation may commence. If the fields perform a random 



wa 
in 



k, the dimensionality of ESLs can be smaller (d = D— 2). Thus, trapped inflation, as envisioned 
^, requires ESLs of dimensionality close to D. 
However, as we shall see in Sec illl Bl a compelling type of trapped inflation occurs naturally if 
D is large and fi ~ x due to the presence of a terminal velocity, even if d = 0. 



B. x ~ /x, a Terminal Velocity on the Landscape 

Consider that the field-trajectory encounters a steep region on the landscape, unsuitable for 
slow-roll inflation. The classical potential leads to an increasing speed so that v = \<p\ becomes 
large. Such an increase in v causes the effective trapping radius /i = y/v/g to grow. For simplicity, 
we consider ESPs only, that is d = so that it is extremely unlikely to run head on into an ESP 
for large D (a generalization to ESLs is straightforward as long as fl > d). 

Once \l ~ x, the trajectory is bound to come within range of neighboring ESPs. For example, 
if D = 1, the inflaton becomes susceptible to 2 ESPs, one ahead and one behind; if D is large, the 
trajectory becomes within reach of many ESPs, of order 2 D . 

Around each of these ESPs, particles are produced (see Sec. lIIIB~T1 for the computation). As the 
inflatons continue to roll, there is an increasing number of particles in the wake of the trajectory 
4 . As we shall derive below, the effect of each individual ESP is small and no actual trapping 
occurs if D is large, but the combined backreaction of all produced particles leads to a resistive 
force opposite to (p. This situation is comparable to that of friction/viscosity in a fluid. 

4 The particles are subsequently diluted due to the expansion of the universe and also back-scatter and decay, see 
Sec. IIIIB1I 
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For large D, this friction increases sharply once \i ~ x, that is around 

v ~ gx 2 , (24) 

preventing a further increase in v and leading to a terminal velocity vt on the landscape, primarily 
determined by the density of ESPs and the coupling strength between the inflaton condensate and 
the new light degrees of freedom. Thus, depending on x and g, inflatons could evolve slowly even if 
the potential is steep. Since no single ESP has a strong influence on <p, we expect the trajectory to 
be largely unchanged, that is, <p/v still follows the potential gradient, but with a maximal speed vt 
determined by balancing the drag force, not Hubble friction, with the driving force of the potential. 
A similar argument can be made for ESLs as long as d <C D. Note that this terminal velocity 
is present at weak coupling, opposite to the speed limit at strong coupling in 
DBI-inflation [33]. 



32] which led to 



We conclude that driving inflation on the landscape may not require a flat potential after all. 
If ESLs are ubiquitous, fields evolve slowly due to quantum backreaction. This type of inflation is 
a generalization of one-dimensional trapped inflation, with the important difference that no actual 
trapping occurs. Furthermore, whereas ESPs need to be aligned neatly like pearls on a string if 

n 

D = 1, which may or may not be natural [18], trapped inflation can result for randomly distributed 
ESPs in higher dimensions. Further, the terminal velocity becomes independent of the slope of the 
potential in the large D limit. 

In the following we analyze the physical processes behind higher dimensional trapped inflation 
in more detail, before providing a phenomenological model and computing some observational 
consequences. 



1. Particle Production near ESPs and Computation of the Terminal Velocity 

Let us consider particle production at a single ESP which the trajectory passes by at a distance 
& t ^esp with constant velocity v (see Sec lIIBI for the notation) 5 . Ignoring backreaction onto 
the trajectory while (f> is close to the ESP, we can compute the number density of x particles with 



wave- number k after the encounter, t > i F ,sp, to 



, k 2 + g 2 fi 2 \ f a(t ESP ) \ 3 

n * =fflq,, ~ 7r— ^J(-<o"J ' (25) 



5 Since particle production near an ESP is exponentially suppressed the larger the distance to the ESP, the majority 
of particles are produced during a short time-frame, when the distance is minimal; this justifies treating v as a 
constant while particles are produced. 
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where we incorporated the subsequent dilution for t > t ESP due to the expansion of the universe. 
The corresponding energy density is 

J^Y n k^/k 2 + 9 2 (<p(t) - ipESp) 2 ~ 9 W(t) - (pesp\ n x , (26) 



where 0esp is the location of the ESP, \<p(t) — (pesp\ = vVo + y2 ~~ *esp) 2 ; an d 

nx = / ** K M!£ e -^/. ( %> )' . (27 ) 



(2tt) 3 re (2vr) 3 V a(t) 

We assume that x-particles are sufficiently stable to ignore their decay, but the unavoidable scat- 
tering off the inflaton condensate lowers the number density over time by an additional factor of 
at most exp(— g 2 n x /(2fj,Qv)) < 1 |8(. For simplicity, we ignore this factor in the following 6 . 

Before we compute the backreaction onto <p, let us investigate up to which distance we need 
to consider ESPs. If we find ourself at a random position in field space, one might guess that 
only nearest neighbors roughly a distance x away need to be considered. However, there are two 
effects which render the effective impact parameter bigger than x. Consider a D-dimensional 
sphere of radius I in field space with an ESP density a = l/x D ; solving Ve>I d <J = 1 gives the radius 
above which we expect at least one ESP inside the sphere. Since V D = ir D / 2 /T{D /2 + 1) and using 
the asymptotic form of the Gamma function for large arguments, T(x) ~ ^/2nx x ~ 1 ^ 2 exp(— x), we 
arrive at 



,SS *V^- (28) 

considerably larger than x for D>1. 

The effect of ESPs that are further away is exponentially suppressed according to (l27|) . but 
their increased number can counterbalance this suppression to some degree. As a consequence, we 
need to consider not only nearest neighbor ESPs. Let us compare the energy density of x-particles 
at ESPs that were produced at t ESP up to a distance of al perpendicular to the trajectory; we are 
interested in their energy density when (t — t ESP )v ^> al, see Fig. HI Noting that for large D almost 
all ESPs are separated from <p(t ESP ) perpendicular to the trajectory and that the majority of ESPs 
can be found near the rim of the cylinder with radius al in Fig. HI we can approximate the ratio 
of the x-P ar ti c l es ' energy density near the rim to the one of particles at nearest neighbor ESPs by 

pfT n i / irql 2 (a 2 - 1)\ 

Ptot _ a D-i exp — y_A L . (29) 



^nearest 
"tot 



The presence of a terminal velocity and even its value in the large D limit does not change if this factor is kept. 
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FIG. 4: The energy density p^ otal is dominated by x-particles in the proximity of ESPs near the boundary 
of a cylinder of radius (1q = od around the classical trajectory <p(t). 

This ratio is maximal for 

a H rp, if „ >2! i, (30) 

^ y gx 2 ' — e ' 

where we ignored terms of order 1/D and used I 2 ~ x 2 D / (2Tte). Thus, the energy density is 
dominated by x-particles a distance 

Ho = al (31) 

away, which sets the characteristic impact parameter for the majority of relevant ESPs. 

Next, we compute the backreaction onto the inflationary trajectory in the Hartree approxima- 



tion [8l4l0l|. We first note that the combined effect of ESPs in the directions perpendicular to the 
trajectory becomes exceedingly suppressed with increasing D: since [iq oc y/~D, the distribution of 
ESPs near the rim of the cylinder in Fig. d] can be well approximated by a continuous ESP den- 
sity; this approximation is similar to treating electric charges distributed on a macroscopic metal 
cylinder as a continuous surface charge density. Then, due to the symmetry of the configuration, 
the net force perpendicular to the trajectory vanishes, but a resistive force opposite to (p remains. 
We assume that the potential is such that the trajectory is straight over times-scales of interest 7 , 
that is over At > 1/(3H) so that (a(t - At)/a(t)) 3 < 1/e. Denoting with ip x the field along the 
trajectory, the equation of motion for <p% becomes 

i p 1 + 3H<p 1 + ^ + g 2 ( X )\i = 0, (33) 



7 If \ip\ — vt during the last sixty e-folds of inflation, the COBE normalization imposes H w 0.018^/W, see forward 
to (I59p . We deduce that the trajectory needs to be straight over Atp > 19^/vt, which should be sub-Planckian. 
Hence, the terminal velocity needs to be 

vt < 10" 4 (32) 

for internal consistency. 
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where 

.D-l-D 



<T(X)>1 ~ / gn x V D ^- v x- u d^ x (34) 

(35) 



M 5/2 1 ( 2MT/ ^ 



(2vr)3 3^ ^ v "»f-o/V X D 
with /io from (f3Tj) and we took f « const as well as -ff w const, so that (a(t ESP ) /a(t)) 3 dt ESP w 
1/3H. For v < there is a strong exponential suppression oc exp(— Dgx 2 /(2ev )), so that back- 
reaction is negligible and trapped inflation is impossible, in agreement with our prior heuristic 
arguments. We recognize 

i> c = — (36) 

as the aforementioned critical velocity of order gx 2 above which backreaction from particles near 
ESPs becomes important. In the regime v > v c the drag increases sharply: assuming 8 \(p\\ <C 
3if|0i| <C \dV/d(px\ in ([33]) . we can compute the terminal velocity by equating the drag with the 
driving force caused by the potential gradient to 

vt = \<pi\& gx 2 A (38) 

where we defined 

Note that A — > 1 in the limit of D —> oo; as a consequence, the terminal velocity becomes 
independent of the slope of the potential. Further, keeping only the leading order terms in 1/D, 
we recognize that the terminal velocity is simply 

v t « gx 2 (40) 

in the large D limit, in agreement with our heuristic arguments of section Ull Bl This is a remark- 
ably simple result, providing a general speed limit on the landscape entirely determined by the 
distribution of ESPs and the coupling strength of the additional light degrees of freedom to the 
moduli, but independent of the slope of the potential. 

If the classical trajectory is intertwined on small scales (Aip < 19y/vt) the above derivation is 
inapplicable; we comment on this regime in Sec. IIII EL 



Using the COBE normalization, see forward to (J59j , yields H ~ 0.018^/vt and the condition 3-ff[y>i| <C \dV/dtpi\ 
imposes an upper bound on the terminal velocity of 

v 3 / 2 < 18.5\dV/d<pi\ , (37) 

indicating that the potential needs to be steep enough for v t to provide a speed limit. In essence, this is the 
opposite of the requirement for a potential to support slow-roll inflation. 
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2. Different Regimes of Trapped Inflation 



Given a reasonably straight classical trajectory, there are two distinct regimes of inflation that 
could occur once v becomes of order u t : 

1. The classical trajectory is followed with a constant velocity close to vt; this requires a large 
D to prevent overshooting and extended oscillations around the terminal velocity. 

2. The classical trajectory is followed with a velocity that oscillates around vt, never settling 
into the constant velocity regime before inflation ends. 

To decide which type of inflation takes place, let us estimate the relevant time scale at which 
backreaction of particles near ESPs becomes stronger than the classical force due to the potential. 
In this section we neglect the expansion of the universe; as the velocity approaches vt from below 
we expand 



v « v t + CL\t 



(41) 



before backreaction becomes important; here, we defined a± = \dV/dipi\ t and set t ESP = 
(particles start to be produced at t = 0). Hence, the characteristic time-scale for the classical force 
is 



v t 



On the other hand, the force due to backreaction increases in time as 



g 2 (x) 2 



Mo) 



gn x V D -i(j,Q x x D dip! 



{g{v t + ai t)fli _ j 
(2vr) 3 6 

(9Vt) 5/2 (. . t\ 
(2vr) 3 V t c J 



2tt x 



D 



(D+4)/2 



D + lt 

D x ' 



(42) 

(43) 
(44) 
(45) 



where we used Vd-i^q 1 /x d a D / fJ,o\/ (D + l)/(27r), a = \J ev/(gx 2 ) from (f30l) and hq from 
9 2 {x) tpi =\dV/8i Pl \ 

*esp an( ^ s °l v i n S f° r t leads to the characteristic time-scale for 



3B- Setting 



backreaction to dominate over the classical force, 



D 



2 ( D + 4 
LambertW 

V 2 

2 (^f 

tESP (^) 6 / 



L> + 4 
2 



dV 



*ESP 



D (2tt) 6 
v t V D + 1 (gvtfl 2 





dV 




dipi 



(46) 
(47) 
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Here, we used (. . .) 2 /(-°+ 4 ) ss 1 + ln(. . -)2/(D + 4) before solving for t b ackr.> and in the last step we 
kept leading order terms in 1/D only and used vt = gx 2 as well as the asymptotic form of the 
LambertW function for large arguments, LambertW(x) ~ ln(x). For 

ibackr. ^ tel. (48) 

backreaction dominates before the classical potential has a chance to significantly alter the velocity. 
As a consequence, the classical trajectory is followed with a constant velocity close to vt once v 
reaches Vt , without prolonged oscillations around vt ■ (|48p along with (|46p imposes a lower bound on 
the dimensionality of field space, D S> 2\n(\dV/ dipi\ 2 ESp (27r) 3 / (gx) 6 ) . If we take a steep potential 
\dV/d<pi\ 

t E sp ~ ^' iar § e ^ couplings g ~ 0.1 and ubiquitous ESPs x ^$ 10 so that vt ^$ 10 and 
([32]) is satisfied, we get D 3> 94. If the potential happens to be shallower, as in Sec lIII C\ lower 
values of D suffice, whereas smaller values of gx require slightly bigger values of D. Note that 
regardless of the slope, we need at least D > 10, since we ignored terms of order 1/D throughout. 

If D is less then the minimal value imposed by (|48p , one might still get inflation of the second 
type: the velocity overshoots, even more light particles are produced since /xo increases, backreac- 
tion sets in and the velocity is forced 10 again well below vt, mimicking the speed of cars during 
rush hour. Due to the changing velocity, we expect bumps in the power-spectrum as well as strong 
additional non-Gaussianities, rendering this regime less appealing. We leave this class of models 
to future studies. 

C. Phenomenology of Trapped Inflation in Higher Dimensions, v ~ const 

Consider a simple multi- field inflationary model, V = ^\ Vi with Vi = m 2 (p 2 /2 and 
fi{pini) = (Pj(Uni) f° r all i,j = 1 . . . D. If we denote the effective field along the trajectory with 
<f = (%2i V 9 ?) 1 ^ 2 ' we § e t V{<p) = iTT' 2t f 2 /2. Let us consider a regime that does not yield slow-roll 
inflation, ipi n i = <p(ti n i) = 1 so that field values remain sub-Planckian (we use 87rG = 1 if not 
stated otherwise in this section). If the trajectory is in a D-dimensional space (D 3> 1) with 
ubiquitous ESPs (or more generally ESLs), we concluded in section fill B II that the speed of ip is 
bounded from above by the terminal velocity vt, see equation (|38p . Furthermore, in the large D 
limit vt Ri gx 2 , independent of the slope of the potential; as a working example that satisfies all 

9 Constraints from IR cascading require g < 0.2, see Sec lIII DTI 

10 The velocity may actually turn around so that the trajectory becomes temporarily trapped near 9?(tesp) up until 
X particles are diluted sufficiently by the expansion of the universe and v increases again. 
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constraints in this paper, we choose g = 0.1 and x = 10~ 3 ' 5 ; then the coupling constant is below 
the bound originating from IR-cascading (see forward to Sec. IIIID ip . vt = 10~ 8 is well below the 
bound in f|32|) originating from the requirement that the distance over which the trajectory has to 
be straight be sub-Planckian, inflation lasts long enough, see forward to (|55p . and (|37p is satisfied 
so that Hubble friction is subdominant. The dimensionality of moduli space needs to be above the 
bound imposed by (|48p to guarantee that ip rolls at the terminal velocity 11 , but the value of D 
does not enter directly into observables as long as it is large enough (D •> 37 is sufficient in this 
section). The Hubble slow evolution parameter becomes 

# 3v t 1 
£ 55 -JP = W^ 6.5x10-^ ' (50) 

where we used the COBE normalization in (159p to set the inflationary scale. Since e should be 



much smaller than one, we find an upper bound for the terminal velocity 

v t < 10" 4 , (51) 

which coincides with the one in (|32p . Inflation takes place up until e ~ 1, that is for <p > ip en d ~ 
V3vt/m ps 0.004 < 1. The energy density at this instant is p = 2v\ = 2 x 10 14 , corresponding to 
an upper bound on the reheating temperature 12 of T r < (30p/ (S^NrG)) 1 ^ ~ 8.8 x 10 11 GeV where 



1PJ and reinstated 



we used N r ~ 1000 as a estimate for the number of relativistic degrees of freedom 
8ttG. To put this another way, to guarantee a high enough temperature for nucleosynthesis, say 
T > lOMeV, we need a terminal velocity of at least 

v t > 10" 35 , (52) 

which is a mild lower bound on vt- 



11 Since the COBE normalization sets the inflationary scale in (|59[) , the mass becomes m m v6 x 0.018-yiit for 
fini ~ 1; then the bound on D imposed by Q48p becomes D S> 37 since dV/dip = m 2 tp ini w 0.0019w t . Note that 
the slope satisfies (dV/ dip) / (3H) ~ ra^J (2/3) > vt so that the terminal velocity is indeed a speed limit. One can 
further check that 

Px « v 2 t /2 « Ptfi , (49) 

as required for inflation to take place (x-particles dominate the equation of motion for ip, but their energy remains 
subdominant.). 

12 The actual reheating temperature depends on several factors: if preheating via parametric resonances (unlikely, 
see Sec lB 2p and/or tachyonic instabilities takes place, the reheating temperature is close to this upper bound, 
which puts it at odds with bounds originating from thermal relics, such as gravitinos (T r ^ (10 s — 10 8 ) GeV, see 
[ifljEtJ or [5ll ] for a brief review). If reheating takes place via the old theory of reheating [IFj, H3, 1 the inflatons 
continue to oscillate around ip — until H becomes comparable to their decay rate(s) V. An estimate of the 
reheating temperature is then given by T r w (90r 2 /(87r 3 iV r G)) 1 . 
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The number of e-folds becomes 



N= l Hdt = I ^^-^fm^i, (53) 

Jt end J Vend V U *2x/6 

where we neglected ip en d <C Vini- As a consequence, we can write 

2 

~ fa« (54) 

Using the inflationary Hubble scale from (|59|) and t/3j n j = 1 in (|53|) . we get the so far strongest 
upper bound on the terminal velocity 

vt < 1(T 8 (55) 

from demanding that N 60 in order to solve the usual problems (flatness, horizon, etc.) of the 
big bang. We chose the inter-ESP distance x = 10~ 3 ' 5 in hindsight to satisfy this observational 
constraint. For different potentials, the numerical pre-factor in (|54p changes, but the scaling 
e oc ipj ni /N 2 remains. Note that e is smaller than in the slow-roll case where esr oc ip 2 ni /N. 



1. Observational Consequences 

Let us focus on adiabatic perturbations, which are properl y d escribed by (the Fourier modes 
of) the Mukhanov variable vt = zC,t with z = aip/H (see 12| for a review). Here £ is the 
curvature perturbation on uniform density surfaces. We are interested in evaluating its power- 
spectrum, which gets imprinted onto fluctuations of the cosmic microwave background radiation 
after inflation. If e evolves slowly, which is the case during inflation, we can approximate 

a(r) oc (-r)-( 1+£ ) , (56) 

where r = — oo ... is conformal time, adr = dt and d( )/dr = ( )', solve the equation of motion 
for v k 

v£+(k 2 -^)v k = (57) 

analytically in terms of Hankel functions Q] , match to the Bunch Davies vacuum in the far 
past (vfe = exp(—ikr)/V2k), expand the solution on large scales and translate back to Cfc- 
In the end, the amplitude of the power-spectrum at horizon crossing (k = aH) becomes 

#2X2 



"<- [ .**L* (58) 
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to leading order in e, which is set by the COBE normalization Pt = (2.41 ± 0.11) x 10 9 |55j. This 
means we need an inflationary scale of 



fl taf .«0.018 v ^. 



(59) 



Further, the scalar spectral index reads 

n, - 1 



d lnV c 
dink 



3- 2i/ 



2Af 2 



(60) 
(61) 



where v is defined by 



z^ i> 

z 



2 _ 1 
4 



(62) 



and we used z"/z ~ 2(1 + 3e)/r 2 . The power-spectrum is closer to scale invariance than in the 
slow-roll case, since e oc (p 2 ni /N 2 . 

The power-spectrum of gravitational waves has an amplitude and spectral index of 1 1 . 

2H 2 



7T- 1 



which leads to the usual expression for the tensor to scalar ratio r 



(63) 
(64) 



r = — - = 16e . 



(65) 



The small va. 
as Planck 



ue of e renders r unobservably small with the current generation of experiments, such 



2. Discussion 

In the simplest realization of trapped inflation in higher dimensions with v = vt = const, the 
Hubble slow evolution parameter is suppressed by an additional power of ./V (the number of e-folds) 
if compared to a slow-roll setup. As a consequence, the scalar power-spectrum is closer to a scale 
invariant one, in tension with the WMAP5 results 55J at the 2a level, and the tensor to scalar 



ratio r is unobservably small. Thus, one could firmly rule out this model of trapped inflation in 
higher dimensions if error bars around n s < 1 tighten or primordial gravitational waves of order 
r ~ 10 -2 be observed. 
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However, we worked with a simple potential resulting in a straight trajectory. In general, 
one might not expect such a smooth potential on the landscape, but one leading to a curved 
trajectory. In such a case we can not neglect isocurvature perturbations, since they feed into the 
adiabatic mode with every twist and turn, leading to additional deviations from scale invariance 
and potentially observable non-Gaussianities. Further, if the trajectory twists and turns on small 
scales, our derivation of the terminal velocity is not applicable any more. We comment on this 
regime in Sec lIII El 

We would like to emphasize that the potential does not influence the speed of the effective field 
along the trajectory in trapped inflation if D is large. Hence, less tuning of the potential is needed, 
but we require a high dimensionality of field space (so that backreaction increases sharply at vt 
and (|15]) is satisfied) and ubiquitous ESPs/ESLs (x < 10 -3 ' 5 , so that vt is sufficiently small, and 
inflation lasts long enough, see ()55j) ) . Further, we impose g < 0.2 to prevent perturbations from 
IR-cascading to dominate, see Sec. IIIID II Other conditions imposed throughout are 

• the trajectory should not be strongly curved over a distance covered in a few Hubble times, 

(EH), 

• the slope of the potential needs to be large enough to guarantee that Hubble friction can be 
ignored, ([37]), 

• the potential energy needs to dominate over the kinetic one, which in turn needs to dominate 
over p x , to guarantee that the universe inflates (|49p . 

• the reheating temperature needs to be large enough for nucleosynthesis to take place, (|52p . 

These requirements are all satisfied in the simple model presented in Sec. IIII CI 

Another advantage of trapped inflation over ordinary slow-roll inflation is its independence on 
the initial velocity, a quantity that requires fine tuning in many slow-roll models. If the initial 
velocity happens to be bigger than vt, a plethora of light states is produced rapidly whose backre- 
action causes a sudden decrease of v well below vt, potentially leading to temporary trapping near 
the initial location. Once x-P ar ticl es are diluted by the expansion of the universe, the speed picks 
up again due to the classical potential and settles into the constant velocity regime of trapped 
inflation if D is large enough. 
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D. Additional Observational Signatures 



We would like to comment on two additional observational signatures brought forth by particle 
production during inflation. 



1. IR-Cascading 



A crucial aspect of trapped inflation is the ongoing production of particles during inflation 
near ESPs. As x-P ar ^ c ^ es become heavy, they re-scatter off the homogeneous inflaton condensate 
and cause a cascade of inflaton fluctuations in the IR 



34, 



35], leading to observable signatures in 



the power-spectrum. For particle production at a single ESP, an additional bump like feature at 



H results, which can be approximated by 



34, 



m , /7re\ 3 / 2 / k V 



(66) 



where the amplitude scales with the coupling g; for example, comparison with lattice field theory 
in a simple, one dimensional m 2 (p 2 inflationary model yields ^4 ESP — f.01 x 10~ 6 g 15//4 for g 2 = 
1,0.1,0.01. Note that the shape in (|66p is only reliable for 10 -7 ~ g 2 ~ 1. Bumps with an 
amplitude of up to 10% of the primordial power-spectrum on scales probed by current CMBR 



35], leading to an upper limit on the coupling 



experiments are still consistent with observations 
of 5 2 <10- 2 . 

However, in the case of trapped inflation as investigated in Sec JIII Cl particles are produced 
at a multitude of ESPs. The resulting superposition of bumps leads to a constant contribution 
to the power-spectrum. It is expected that this contribution, albeit scale invariant, is strongly 



non-Gaussian and thus ruled out as the dominant one 
In 



35] 



13 



35] observational bounds on the amplitude of features were given under the assumption 
that all ESPs lie on the trajectory, as in one dimensional monodromy inflation [18], resulting in 
a superposition of bumps with identical amplitudes. For a large ESP or bump density (hi = 
kxe^ A with A < 1), it was found that the constant contribution scales as A ; consequently, 
the maximally allowed amplitude is proportional to A, roughly A < 2A x 10 -9 , which provides 
and upper bound on g 2 . 



13 A computation of Non-Gaussianities from IR cascading is in preparation by the authors of 3S]; observational 
bounds on non-Gaussianities might yield stronger constraints on particle production during inflation than consid- 
erations based on the power-spectrum. 
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This may sound alarming, since trapped inflation in higher dimensions relies on a multitude 
of ESPs, roughly n ESP = Vd/Zq 3 /x d ~ a D ~ e D / 2 3> 1, at any given time. There is however a 
crucial difference, rendering the conclusions of [33] only indirectly applicable. The x-P ar ticl es that 



dominate the energy density p x in Sec HUB II are produced at ESPs a distance po = \/D/(2ir)x 
away from the trajectory. As a consequence, particle production at any single ESP is weak; to be 
concrete, the particle number at each ESP is exponentially suppressed 

( k'+g^l 
rik = exp — 7r 



ex P ( -^2) exp (-§ ) ■ (67) 



The amplitude of each bump should carry the same exponential suppression oc e~ D l 2 . Hence, 
the effect of n ESP ~ e D l 2 ESPs at any given time leads to a constant contribution to the power- 
spectrum with a magnitude corresponding to the amplitude of a single feature, A ~ 10~ 6 g 15 / 4 ; 
then the upper bound of order A < 2 x 10 -9 yields the relatively weak bound of g % 0.2. Of 
course, this heuristic argument should be followed up with a more rigorous computation, which we 
leave to a future publication. 

We would like to point out that large values of g are not required for trapped inflation in higher 
dimensions to operate; in fact, lower values of g might be desirable, since smaller terminal velocities 
can be achieved with less ubiquitous ESPs. 



2. Individual Massive x-particles 

Throughout trapped inflation in higher dimensions, X"P ar tid es with a (^-dependent mass are 
present. The larger the mass of a particle, the further away from the current position on the 
trajectory is the ESP at which the particle has been produced. As this distance in field space grows, 
particles at a given ESP become more and more diluted due to inflation while their mass increases. 
As particles become more massive, their decay rate to lighter degrees of freedom increases, T x oc m 
(i.e. T x = h?m/{&i:) in the presence of a coupling to a spinor field —hipipx [ill)) an d they decay 
once T x becomes of order H. If they decay while their density is still large enough so that they 
can be described by a homogeneous field (as we assumed throught the paper), they will leave no 
additional observational traces, since p x <C <C p v . Note that they need to be stable for at 
least a few Hubble times after their creation in order for backreaction to become important in the 
first place. 
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On the other hand, if ^-particles are stable enough that their number density becomes exceed- 
ingly low before they decay, we have to consider the effect of individual, massive particles (a few 
per Hubble-horizon) that are directly coupled to the inflatons. In the vicinity (in real space) of 
such a particle the inflatons are held back. As a result, the energy-density in a spherical re- 
gion around the particle is slightly increased. This in turn delays reheating and causes a circular 
cold spot in the CMBR with an amplitude that is proportional to the induced over-density 14 , see 



i.e. 



36-38]. Since these massive x-P ar ticl es are resupplied throughout trapped inflation, we arrive 
at an additional source of perturbations, consisting of the superpositions of spherical over-densities. 
We expect these additional signatures to be ill suited to provide the dominant contribution to the 
power-spectrum. An investigation of this interesting effect is in preparation. 



E. A Bending Field Trajectory 

The force due to backreaction of x-P ar ticles is not anti-parallel to the velocity if the trajectory 
is curved on length-scales of order vt/(3H), where 1/(3H) is the time-scale over which backreaction 
of particles at a given ESP is important. As a result, any bending of the trajectory is enhanced, 
causing the trajectory to twist and turn. Inflation may still take place, but since the trajectory 
bends considerably during any given Hubble time, isocurvature fluctuations constantly feed into the 
adiabatic mode, making the power-spectrum less scale invariant and enhancing non-Gaussianities. 
In order to say anything more concrete, we need to specify the classical potential that determines 
the classical trajectory. 



1. A Meandering Path 

Assume that the inflaton potential is smooth on large scales with a constant average gradient in 
the direction of tpi, but with random small scale bumps and dips on top. Whenever the trajectory 
encounters a bump it gets disturbed, causing it to find a way around the bump. The resulting 
tra ject ory is reminiscent of the path a meandering river takes flowing down a mountain side, see 



i.e 



57 



591 ] for related proposals. 
On such a potential, the trajectory encounters many bifurcation junctions: for example, if 
D = 2, the trajectory could go to the left or right of a bump. Which way is chosen depends 
sensitively on the initial conditions. If we deal with many bumps, the displacement perpendicular 



We would like to thank A. Dahlen for pointing out this possibility. 
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to ip i can be described by a random walk, whereas ip\ remains dominated by the averaged gradient. 
If D > 4 the probability that trajectories reunite approaches zero, since a random walk is not self 
intersecting in three or more dimensions. 

Inflation in such a model is troublesome: at best, the universe is covered by domain walls after 
inflation if by some miracle in all the different inflationary Hubble patches that are observable 
today inflation terminated and reheating commenced simultaneously to within a fraction of a 
Hubble time, even though the trajectories in each patch are different. Without such a miracle, 
inflation ends at different times, leading to pronounced hot/cold spots in the CMB, that are not 
observed. 

Adding ESLs and backreaction to such a model enhances these problems. Therefore, we do not 
consider meandering trajectories further. 



2. A Curved Path 



Due to the problems caused by bifurcation junctions, such as domain walls and pronounced 
hot/cold spots in the CMBR, we would like to shift attention to potentials that channel the tra- 
jectory along a single path throughout the last sixty e-folds of inflation. This means that we would 
like to investigate potentials that have a positive semideflnite slope in all directions perpendicular 
to the trajectory, see for example 



60] for a recent proposal of an in-spiraling trajectory in two 
dimensions. 

If we consider models with a clear channel, that is with slopes within the channel walls that 
are bigger than the slope along the channel, and if the classical trajectory is only curved but 
not intertwined over length-scales of interest, we expect the effect of backreaction to be still most 
pronounced in the direction opposite to the velocity. Thus, we still expect the presence of a terminal 
velocity, as derived in Sec. IIIIB 11 If backreaction alters the trajectory drastically, i.e. by forcing 
the fields over a channel- wall, viable inflationary models could still result, but little can be said 
generically. We leave a numerical investigation of concrete models to future studies. 



IV. CONCLUSION 



Whether or not moduli get trapped at extra species loci (ESL), which are often associated with 
additional symmetries, has important consequences for selecting vacua on the landscape. If the 
trajectory in field space comes close to such a locus, particles are produced and backreaction can 



29 



trap the trajectory. In this paper we gave a simple geometric argument to highlight that this type 
of moduli trapping is unlikely if the dimensionality of moduli space is considerably larger than 
the dimensionality of the loci; for D > d + 1 (or D > d + 2 if the trajectory is intertwined), the 
probability of getting trapped is suppressed by (fi/x) D ~ d ~ 1 (or (fi/x) D ~ d ~ 2 if the trajectory is 
intertwined), where x the characteristic inter-ESL distance, fi = y 'v/ 1 g < x is an effective trapping 
radius, v the velocity in field space, and g a coupling constant. We concluded that the sole presence 
of ESLs does not guarantee a dynamical preference of high symmetry vacua on the landscape, but 
since the effective trapping radius increases with the velocity, ESLs are always crucial at large 
speeds, for instance if moduli are responsible for inflation. 

For D » 1 and d = (extra species points or ESPs, a generalization to ESLs is straightforward), 
we derived a terminal velocity vt = g 2 x on the landscape: once the speed approaches vt, the 
trajectory comes within reach of a plethora of ESPs at which x-particles are produced. Their 
combined backreaction acts opposite to the velocity if the classical trajectory is straight over 
distances covered in a few Hubble times (we allow for a classical potential of the moduli). Since 
backreaction increases sharply for v > vt if the dimensionality of moduli space is large, vt acts as 
a speed limit. The classical trajectory is followed with this speed as long as the slope along the 
classical potential is large and the path is reasonably straight. 

The presence of a terminal velocity offers the intriguing possibility to drive inflation with po- 
tentials that do not otherwise support slow-roll inflation. This type of inflation is a generalization 
of trapped inflation. The main difference is the independence of the speed from the slope of the 
potential, as long as D is large enough. 

We then provided a phenomenological model of trapped inflation in higher dimensions, assuming 
a simple quadratic potential and ubiquitous ESPs (x < 10~ 3 5 , g ~ 0.1). For these parameters, more 
than sixty e-folds for inflation result for sub-Planckian excursions in field space and observations 
are within observational bounds (the power-spectrum of curvature perturbations is slightly red 
and gravitational waves are suppressed). We would like to highlight that no fine tuning of initial 
conditions is needed and constraints on the potential are lessened (the COBE normalization still 
sets the inflationary scale, but the 77-problem is alleviated). Additional observational signatures 
are expected, such as an extra contribution to the power-spectrum from IR-cascading (less than 
a few percent if g < 0.2), potentially large non-Gaussianities and additional circular cold-spots in 
the CMBR caused by individual x-particles that became heavy before being diluted away. We plan 
to come back to these interesting signatures in a future study. 

In this paper, we took a bottom-up approach and parametrized the undoubtedly rich structure 
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of the landscape and the distribution of ESLs by a few characteristic parameters, such as x, that 
we left undetermined. Concrete implementations within string theory are needed to see whether 
or not the distribution of ESLs and the potential on the landscape supports trapped inflation in 
higher dimensions. It appears that ubiquitous ESLs and steep potentials are common, but we leave 
an in depth analysis to future studies. 
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Appendix A: Trapping Probability, Improved Estimate 

It is instructive to estimate the trapping probability by means of a volume argument, which 
does not require any artificial shifting of ESLs. Consider the volume in the vicinity of the classical 
trajectory such that the shortest distance to (p(t) is less than \x 

V d = sV D ^ D - 1 . (Al) 

If we consider ESPs (d=0), we can estimate the probability of trapping simply as p c i ~ V c i/x D = 
Vd^i(s / x){[i / x) ^ 1 . Note that the numerical pre-factor differs by 1/D if we compare the prob- 
ability with ([20]); the origin of 1/D in ([20]) can be traced to shifting ESPs to the surface, so its 
absence is expected. 

We can generalize this argument again to higher dimensional ESLs, but we need to specify the 
orientation of ESLs in order to derive analytic expressions for the trapping probability 15 : consider 
a periodic grid of ESLs (d > 1) with a grid-spacing of x, such that one extension of an ESL is in 
the same direction as (p{t). For example, if d = 1 and D = 3, we consider a cubic lattice where any 
given extra species lines intersect with two other extra species lines at right angles after a distance 

1 We did not need to specify the orientation in section III C 21 since we treated ESLs as if they were laid out over the 
surface of the hyper-sphere with radius s around 0(0); consequently, section III C 21 is more general with respect to 
the orientation of the ESLs, but our arguments depend on the artificial movement of ESLs to the surface, which 
leads to a lower bound of the probability. 
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ESL 



FIG. 5: Schematic in D = 3, showing V c i = sVd~iH D1 (the gray cylinder), two (exemplary) extra species 
lines (d = 1) separated a distance x that are perpendicular to the classical trajectory (p{t) (not shown) of 
length s; (pit) lies in the middle of the dark gray rectangle (area A± = sn D ~ 1 ~ d VD-i-d = 2/zs), connecting 
the two faces of the cylinder. The probability of getting trapped is p c i rj N±A±/x D ~ d = A^is/x 2 , in 
accordance with (1221). 



x, and one class of lines is parallel to <p(t). Such an arrangement is artificial, but it allows for a 
simple estimate of p c i that provides an approximation for more random distributions of ESLs that 
share the same inter-ESL distance x. 
Given such a grid, there are 

"U-\d-lJ- (d-l)\(D-d)\ [AZ) 
possible orientations of an ESL with one direction parallel to ip(t) and 

possible orientations perpendicular to <p(t). The probability for a parallel (perpendicular) orienta- 
tion to intersect V c i is 

P\\ = Wn-^ D - d Jz- d , (A4) 

1 

x 1 

For example, if d = 1 and D = 3 (so that iV|| = 1 and iVj_ = 2), we can identify 1/x 2 as the 
number density of intersection points between the extra species lines and projections of V c \ onto 



p± = N±s v D _ 1 _ dfJ , D - 1 - d -±- 3 . (A5) 
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planes parallel or perpendicular to <p(t); to be concrete, sVo—x-d^ = Vifis = 2/j.s is the cross 

section of V c i along <p(t); thus, 2fis/x 2 counts the number of intersection points of the ESLs with 
the cross-section. This is identical with the probability of being trapped at one of these ESLs if 
2/j,s/x 2 < 1. See Fig. [5]for a schematic. 

The total trapping probability is given by the sum p c i = p± + P\\, but since \i <C s we can 
approximate 

Pd ~ Px = N ± V D -i-d- (-) ■ (A6) 

Note that this expression is only applicable for p c i < 1. If we compare this probability with the 
lower bound in (|2ip . we observe again the absence of the 1/D suppression (as expected), and a 
numerical pre- factor iVj_, which is determined by the chosen orientation of ESLs along a grid. For 
comparison, if all ESLs were parallel to each other, while having no extent in the direction of (pit), 
the pre-factor would be N± = 1. Ignoring the pre-factor, which needs to be determined from the 
given distribution of ESLs in concrete cases, we find the general scaling 

S / LL \ D — l—d 

Pciocyb-i-d-(-) • (A7) 

x \xJ 

For large D — 1 — d, we observe a strong suppression due to ([j>/x) D ~ 1 ~ d as well as Vo-\-d- We 
conclude that moduli trapping by ESLs is only efficient if d = D — 1, that is ESPs can trap moduli 
efficiently on a line (D = 1), but not on a plane (D = 2), or higher dimensional spaces. 



Appendix B: Examples 



Size Moduli 



It has been argued [7j, [LJ] (see [6l|] for a review) that the size of internal dimensions can be 
stabilized by string gases that contain enhanced symmetry states 16 . String Gas Cosmology (SGC) 
is a low energy field theoretical approach incorporating strings and branes in a cosmological setting 
by means of a gas approximation. If internal dimension come close to the self dual radius, certain 
combinations of winding and momentum modes become massless. Once they are produced in 
sufficient quantities, their backreaction can stabilize internal dimensions temporarily, up until the 
expansion of the universe dilutes them away. 



See [151 ] for a related work in M-theory. 
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There are six internal dimensions in string theory, D = 6, and the light states appear whenever 
any of the six internal dimensions comes close to the self dual radius, d = D — 1 = 5; because 
d = D — 1, this type of trapping is a simple example of an unsuppressed trapping probability. As 
a consequence, temporary trapping appears unavoidable if an internal dimension approaches the 
self dual radius. 



2. Preheating 

Preheating refers to a possible phase of explosive particle production after inflation, whereby a 
large fraction of the energy in the inflaton fields is transferred to other degrees of freedom. There 
are several phenomenological models of preheating employing resonances P, [h| or instabilities 
62, |63] , see for reviews and comprehensive referencing. 

Simple parametric resonance models operate via the mechanism described in Sec. MB} with ipi 
(i = 1 . . . D) identified as the inflatons. These are coupled to a "matter" -field x> which is assumed 
to decay further to ordinary standard model particles. Let's assume that all inflaton fields couple 
to one matter field only, which becomes light at tpi = 0. Because the condition for resonant particle 



i y 

we conclude that all inflatons need to get close to the origin simultaneously (in phase), in order for 
parametric resonance to take place. Even if the inflaton potential has a global minimum at (p = 0, 
it is unlikely that all fields start out in phase and oscillate with the same frequency. Because of 



production is ([15 



the generic dephasing of fields, we expect parametric resonance to be negligib 



couple to the same matter field. This is indeed the case, as shown in 
the example of AA-flation, and by means of lattice simulations in 
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e if several inflatons 



651 ] analytically at 



661 ]. where it was also shown 



that for as little as D = 3, parametric resonance is heavily suppressed (see also 



67]). This is in 



agreement with our conclusion that trapping events are unlikely if d < D — 1, see eqn. ()2ip . One 



can still observe resonances for D = 2 [66j even though = d < D — 2 = 1, since the classical 
potential of the inflatons, which we neglected in deriving (f2"Tj) . leads to repeated attempts of particle 
production at the same ESP 17 . Note that the absence of parametric resonance may be desirable: 



The case D = 2 is even more subtle: in [d9| . [70| it was shown that almost all Fourier modes of % get amplified 
if the oscillating frequencies of the inflatons are not simple ratios. This enhancement is sometime referred to as 
Cantor preheating and was shown to be inefficient for D > 3 in [6q| . see also [67]]. 
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since the reheating temperature is lowered considerably, the overproduction of thermal relics, such 



as gravitinos [4814501] , is avoided. 

In [6a] it was also shown that tachyonic instabilities, that is explosive particle production caused 
by Yukawa interactions of the inflatons to the matter field (£i n t = 'Yli^i l PiX 1 : /2), are left intact if 
the number of fields is increased. Such interactions are not present in the cases we consider in this 
paper. 

Appendix C: Drifting through Moduli Space during Inflation 

During inflation any classical movement of moduli is quickly damped by Hubble friction, as 
we have seen in section III Al However, quantum mechanical fluctuations cause the fields to drift 
in random directions 18 . Consider a single modulus located initially at the origin <p(0) = 0, with 
negligible classical potential (the moduli we consider here are not the degrees of freedom that drive 
inflation; see section ITTT1 for the latter case). Inside the horizon, quantum fluctuations displace the 
field from ip = all the time; these fluctuations become classical shifts of the field after horizon 
crossing. Thus, if we look within our Hubble patch and ask what is the value of <p, homogeneous 



within our Hubble patch, we generically find a non-zero value. If we were to look at an ensemb 
Hubble patches, which inflation creates, we expect to see a Gaussian distribution. Following 



e of 



the probabi 



equation 



68] 



ity distribution of finding a coordinate volume filled with <p ^ satisfies a diffusion 



where T> = H^/Stt 2 , which can easily be solved for H ~ const, yielding the expected Gaussian 
probability distribution 

Pfe^expf-^V (02) 



HH r V H H 
Thus, the expectation value of <p 2 in a Hubble volume is 



(ip 2 ) = J ip 2 P(<p,t)dp (C3) 

= J?*- («) 



18 Here and in the following we assume that Hubble induced corrections to the moduli's mass are negligible. If they 
are present, the fields do not drift further than Sipi ~ H. 
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Since the number of e-folds during inflation is 



N = 



I 



Hdt Ht 



(C5) 



we get the rms value of the modulus after N e-folds 



rms 



H 



(C6) 



If we are not interested in resolving the field trajectory on smaller scales than H/(2n), we can 
approximate the diffusive, continuous drifting of the field in any given Hubble patch by a random 
walk, whereby the field makes a step to the left or right of length H/(2tv) after each e-fold of 
inflation. 

In a higher dimensional moduli-space, each of the fields diffuses independently of the other ones, 
causing a random walk in the higher dimensional space; after each e-fold, a step of length I = H/2ir 
is performed in a random direction. To keep things simple, we approximate the continuous choice 
of direction by a discrete one along a D-dimensional cubic lattice along the (/^-directions. This 
approximation is viable if we do not resolve the path on scales comparable or smaller than the step 
size. 

Given such a random walk in field space 19 , there are several effects determining the trapping 
probability: 

• The maximal length of the random walk in D-dimensions, which is given by the largest 
radius of gyration R\ (see section ICTj) . Using s = s rw = 2R\ in ([22]) yields a lower bound 
on the trapping probability. 

• The shape and volume coverage of a random walk. Even though the ensemble average over 
many random walks has spherical symmetry in moduli space, an individual random walk is 
far from spherical. In the large D limit, one can derive analytic expressions for all radii of 
gyration (see section IC lj) . Within the volume spanned by the radii of gyration, the random 
walk has an increased probability of getting close to an ESL (the Hausdorff dimension of a 
random walk for D > 2 is two; thus, it is more volume filling than a straight trajectory with 
the same s). 

19 We work with a non-compact field space in this section. If moduli space is compact and the trajectory is ergodic 
[Tlj . one may still use the arguments in this section by treating the field space as a covering space of the moduli 
space. 
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• Expansion due to inflation. In an inflating universe, one might consider not only the history 
of our own Hubble volume, but the ensemble of all Hubble patches that became causally 
disconnected during inflation. Since the moduli perform an independent random walk in each 
causally disconnected Hubble patch, one might be tempted to include an additional volume 
factor that enhances the probability. We argue in section IC 41 that this line of reasoning is 
flawed. 

In the following, we assume y/~Nl > x > /x. 



1. Random Walks in D Dimensions 



Consider a random walk of step-length I = H/(2tt) on a lattice in a D-dimensional space. Given 
a walk consisting of N steps (N is identical to the number of e-folds in our case) we can write the 
position after the k'th step as <p(k), that is the i'th modulus field has a magnitude of ipi(k) after 
the k'th step. A walk's shape can be characterized by the radius of gyration tensor, defined as 

HQ 

= n £ &i(k) - M*0 - (<Pj)rJ > ( C7 ) 

k=l 

in analogy to the moment of inertia tensor in classical mechanics. Here, (.) rw is the average over 
a particular random walk, that is 

1 N 

k=l 

The eigenvalues of (|C7j) . Rj, provide an estimate of the walk's extent in the direction of the 
corresponding eigenvector. One can derive analytic expressions for the ensemble average {-) en in a 
l/D-expansion (which becomes exact in the limit of D — > oo), leading to the radii of gyration 

,2\ ^,2^ + 1 A • 3 



W>e»" l l?fr{ 1 + W) ' (C9) 

up to first order in 1/D (j = 1 . . . D). The associated variance is 

which becomes smaller with increasing D. It should be noted that the error of the 1/D expansion 
is smallest for low j, and radii are better approximated than variances [73;]. However, even for 
D = 3 and keeping terms up to order 1/D, the approximation is accurate enough for our purposes; 

37 




FIG. 6: Visualization of the volume sampled by a random walk in D = 3: the random walk covers an 
area that can be estimated by the two largest radii of gyration A rw = s rw .is rw .2 ~ 4i?ii?2, where s rw ^ are 
given by twice the ensemble average of the corresponding radius of gyration. Moduli get trapped if an ESL 
intersects the volume V rw = A rw VD-2^ D ~ 2 = 2A rw /j, (the grey monolith in the picture). Also shown are 
ESPs (d=0) with a characteristic separation of x, so that p rw ~ V rw /x D — 2s rW: is rw ^fi/x 3 (only applicable 
for p rw < 1). 

mean values of R\ from simulations agree with the ensemble average of the eigenvalue to within 
1%, R\ to within 10% and R\ to within 25%. These errors are consistent with expected corrections 
of order 1/D 2 . 

As evident from the radii of gyration in (|C9p . the shape of an individual random walk is far 
from spherical, since a sphere has Rj = R for all j. One can quantify the amount of asphericity by 



a number that is between zero (a sphere) and one (a stick). For large D, the asphericity approaches 
2/5, indicating that random walks in higher dimensions are neither stick-like, as the classical 
trajectory, nor spherical, as an average over many random walks. Further, they all "look" alike, 
regardless of the dimensionality D. 
The random walk is spread over 




(Cll) 



4 + 2D 
4 + 5D 



(C12) 




(C13) 




1/2 



(C14) 
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Hence, the largest extent of the random walk is given by 

y/N + 1 / 3 \'/ 2 

Srw = ~ 21 ( 1+ 4£)j ■ ( C15 ) 

A naive guess for the volume sampled by a random walk mig ht be Uf =1 ( 

s rw,j)- However, this 

does not properly take into account the volume coverage of a random walk; its fractal dimension, 
that is its Hausdorff dimension, is 2 for all D > 2. Thus, a random walk is not volume filling, 
but it covers an analog of a two dimensional surface with an area of order s rWt \s rWt 2- Since we are 
interested in cases where the random walk comes to within a distance of order \i to an ESL, we 
may define an effective volume of (see Fig. [6]) 

V rw = VD-2V D ~ 2 S r w,lSrw,2 (C16) 

. iV -V D _ 2 6( 1 + A)&til. (C17) 
V rw can be used in estimating the trapping probability, see section IC 31 



2. Trapping a Random Walk 

Moduli trapping, and the magnitude of /i, differ from the process described in section III Bl since 
the movement of the random walk is not driven by kinetic energy, but by quantum fluctuations. 
The description in section IIIBI is based on modeling the trapping event in field theory by the 
inclusion of a classical field x, which acquires a non-zero VEV due to explosive particle production 
once the trajectory in moduli space comes closer than \x to an ESL. In this process, part of the 
kinetic energy in the moduli is transferred to %, which then backreacts onto the <fi. However, if 
the moduli simply drift through moduli space, driven by quantum mechanical fluctuations, this 
process can not be operational. 

Nevertheless, there is still a mechanism that can trap moduli: let's denote the distance from 
an ESP after which a field theoretical description with an additional scalar field x becomes viable 
by fj,ff. Once the trajectory comes closer than fj,j t to the ESP, we also have to consider the effect 
of quantum mechanical fluctuations on x, which increases Xrms = \/< X 2 >■ in addition to this 
quantum mechanical drift, the field is affected by the classical potential 

^ = T* 2 X> 2 - ( C18 ) 

Since ^ ipj « /ij t , a classical force driving x back to the origin results. To be concrete, if the 
usual slow-roll approximation for x is used, which is justified since we start out at x = 0, one finds 
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a modified diffusion equation for the probability distribution of finding x at a nonzero 
VEV, 

dP x _ H 3 d 2 P x 1 d ( dV\ 

dt s^ 2 d x 2 3Hd x V x dxJ ' 1 j 

where the last term originates from the classical force on x- This leads to an equation of motion 
for (x 2 ) 

dtf) _ H* 2 /8V \ _ 



dt Air 2 3H\d X 
Since dV/dx ~ <? 2 ^/tX) we can integrate the above to 



After iV > 3H 2 /(2fi 2 ft ) e-folds of inflation, the value of 

Xrms settles into 



*~-Vs??!$- (C22) 

This in turn introduces a classical force onto the moduli oc dV/dtpi ~ ^Xrms^i^ driving them 
towards the position of the ESP. Once they get closer to the ESP, x can drift out further, increasing 
the classical force on the moduli, etc. In the end, the moduli drift around the ESP, while x 
fluctuates around zero, and none of them can escape. Thus, the moduli are trapped near the ESP 
(a generalization to higher dimensional ESLs is straightforward). 

Interestingly, the effective trapping radius is given by /u = jj,f t , which is necessarily bigger than 
the value in section III Bt since the arguments in section III Bl are also based on a field theoretical 
description. Thus, trapping by ESLs is more efficient for drifting moduli than for classically moving 
ones. If Hf t is bigger than the expected inter-ESL distance x, we conclude that moduli are preferably 
located near ESLs after inflation, given that inflation lasts sufficiently long so that s rw S> fift- 

The value of /x/j is model dependent; for example, if the moduli correspond to distances between 
D-branes and the light extra degrees of freedom correspond to strings stretching between them, 
then the strings are only light, and can thus be modelled in field theory by x, if the branes come 
to within the fundamental string length. This in turn determines 

In the following, we leave \x = \Xft general and take s rw > x > \x to estimate the trapping 
probability if /j,f t turns out to be less than x. 
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3. Trapping Probability 



Using the maximal extent of the random walk s rw from (|C15p in (|22|) yields a lower bound on 
the trapping probability 



Prw > Pd(s rw ) « iVj.Vb-1-d— (-) (C23) 



vr \ 4.D y x \x< 

where I = H/ (2tt) and N = J Hdt ~ i^t is the number of e-folds, identical with the number of 
steps in the random walk. 

Incorporating the increased volume coverage of the random walk leads to a better estimate 
of the trapping probability. Analogously to appendix [Aj we distribute the ESLs along a grid. 
The total probability is now a sum of four parts, p±±_ (all ESL directions are perpendicular to 
s rWi i and s rw> 2), p±u and p\\± (one ESL direction is parallel to s rw> i or s rw ^, while all other ones 
are perpendicular) and (one direction is parallel to s rw ,i and an other one to s rWt 2). Each 
probability is given by the appropriate cross-section of V rw (such that the ESLs appear as points 
in this cross-section) multiplied with l/x D ~ d and the number of possible orientations of an ESL 
within this class. Since \x <C s rWt i, s rw ^, we can approximate the total probability by p rw « pj_ j_. 
If ESLs are distributed along a grid, there are 

possible orientations perpendicular to both s rWt \ and s riUj 2, and the probability becomes 

Prw ~ A r _|__|_Vci_ 2 _ f f/i jD 2 d s rWt is rw> 2 D _ d (C26) 

As before, a different distribution of ESLs would change the pre-factor, but not the general scaling 
with fi/x. A comparison with (|C24|) reveals an enhancement of y/N + 11/ /x, caused by the higher 
volume coverage of the random walk. However, for large D — 2 — d, we still observe a strong 
suppression. We conclude, that moduli trapping at ESLs in D dimensions is expected to be efficient 
during inflation only if d > D — 2; for instance, lines can trap efficiently in three dimensions, while 
points can't. 

The expression for p rw in (|C27p is also valid for a classical, intertwined trajectory, that can be 
approximated by a random walk. In this case \i = y/v/g, N is still the number of steps (but not 
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the number of e-folds) and I the step length, all of which need to be determined from the classical 
trajectory. 



4. Expansion Effects 

If a random walk in moduli space begins in some Hubble patch during inflation, e 3 identical 
copies of this Hubble patch appear after the first step (that is after one e-fold), and the moduli 
perform independent random walks in each patch from then on. 

Consider now a random walk of length s rWt i = x. The volume sampled by a single walk of this 
length is (fCTTl) 

2 

V rw = y ii D ~ 2 V D - 2 ■ (C28) 



Thus, it would take about 



n^= 2 f_;y (C29) 



V rw fL D ~W D 



-2 



independent random walks to cover the entire volume of the D-dimensional sphere of radius x. As 
a consequence, the probability of finding at least one Hubble patch where the moduli get trapped 
at an ESL approaches one, if inflation lasted N a dd > hi(n)/3 e-folds in addition to the number of 
e-folds needed to achieve s rw = x. 

Thus, one might be tempted to conclude that inflation enhances the chance of being trapped 
by ESLs. However, this line of thought is flawed: the volume enhancement factor does not prefer 
ESLs; inflation allows for the population of all vacua that are within reach of a random walk, 
not only those close to an ESL. Hence, the presence of ESLs is irrelevant if arguments to select a 
vacuum are contingent on this indiscriminatory population of moduli space. 

An additional problem with choosing vacua randomly as above is the presence of domain walls 
at late times: if the fields drift through moduli space during the last sixty e-folds of inflation, 
different regions of the universe today have distinct values for the moduli- fields. In this section, 
we ignored any classical potential for the (pi, an approximation that may be justified during high 
scale inflation but certainly not today. As a consequence, our universe would be filled with domain 
walls separating regions with different values for the a setup that is ruled out observationally. 
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